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SUMMARY 


This  document  is  a  final  report  on  the  measurements  of  the  light-scattering  properties  and 
aerodynamic  drag  characteristics  of  blow-off  particles  generated  in  simulated  thermonuclear  blasts. 
The  "thermal  dust  particles"  examined  were  generated  in  experiments  performed  by  Science 
Applications  Inc.  (SAI)  under  a  separate  DNA  contract,  using  a  solar  furnace  in  Odeillo,  France.  The 
samples  obtained  from  SAI  consisted  of  filters,  wail  wipings  and  loose  dust  taken  from  various 
locations  in  the  solar  furnace,  and  they  correspond  to  various  soil  samples  placed  in  the  furnace  for 
the  tests. 

The  data  contained  herein  include  particle  size  and  shape,  complex  refractive  indices,  single 
particle  albedo  results  and  aerodynamic  drag  forces  as  a  function  of  gas  velocity.  A  new  technique 
was  developed  and  applied  for  the  first  time  here  for  the  direct  measurement  of  aerodynamic  drag 
forces  on  micron-size  particles. 

The  optical  properties  of  blow-off  particulate  matter  in  the  aerosol  size  range  were  found  to  be 
similar  to  the  bulk  properties  of  mixed  silicates  with  the  measured  real  components  of  the  complex 
refractive  indices  in  the  range  1.391  to  1.647.  The  wavelengths  used  in  the  three  phases  of  the 
research  were  in  the  range  454.5  nm  to  632.8  nm.  In  this  range  the  refractive  indices  were  so  nearly 
constant  that  no  dispersion  Nation  could  be  calculated  from  the  data  to  correlate  the  refractive  index 
as  a  function  of  wavelength.  The  imaginary  components  of  the  complex  refractive  indices  were  found 
to  be  near  zero  within  the  limits  of  experimental  error,  indicating  that  little  absorption  of  the  light 
occurred  in  the  visible  region  of  the  spectrum  of  electromagnetic  radiation.  Single  particle  albedos 
were  determined  to  be  approximately  unity,  ranging  from  0.68  to  0.99  based  on  the  measured 
complex  refractive  indices. 

Although  many  particles  had  irregular  shapes,  the  drag  forces  on  them  were  approximately  the 
same  as  for  spheres.  With  few  exceptions  the  aerodynamic  drag  characteristics  of  the  particles  were 
found  to  satisfy  the  Stokes  relation  for  the  drag  force,  using  the  aerodynamic  radius  as  the  equivalent 
radius  of  the  sphere. 


PREFACE 


The  work  reported  here  represents  the  third  and  last  phase  of  baseline  studies  on  the  optical 
properties  and  aerodynamic  drag  characteristics  of  blow-off  particles.  The  results  of  Phase  I,  which 
involved  the  development  of  instrumentation,  measurement  techniques  and  computer  software, 
were  summarized  in  Report  No.  DNA  6274F,  dated  18  April  1983.  Phase  I  was  supported  under  DNA 
Contract  No.  DNA  001-80-C-0213.  The  work  on  Phase  I  was  performed  at  the  Institute  of  Paper 
Chemistry  in  Appleton,  Wisconsin  under  a  subcontract  with  the  University  of  New  Mexico.  Phase  II, 
which  involved  light-scattering  measurements,  was  summarized  in  a  report  on  DNA  Contract  No.  DNA 
001-84-C-0224  dated  29  February  1984.  Part  of  the  work  on  Phase  II  was  performed  at  the  Institute 
of  Paper  Chemistry,  and  the  work  was  completed  at  the  University  of  Washington  under  a  subcontract 
with  the  Institute  of  Paper  Chemistry. 
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SECTION  1 
OBJECTIVES 


The  objectives  of  this  research  were  to  measure  the  aerodynamic  drag  characteristics  and 
light-scattering  properties  of  dust  particles  generated  thermally  in  a  solar  furnace.  The  following 
properties  were  measured  and  are  reported  here: 

1 .  The  aerodynamic  drag  force  as  a  function  of  Reynolds  number. 

2.  The  aerodynamic  diameter. 

3.  The  particle  mobility  (ratio  of  velocity  to  drag  force) 

4.  The  light-scattering  phase  function  (scattering  intensity  versus  angle)  as 
a  function  of  wavelength  of  the  incident  radiation. 

5.  The  complex  refractive  index. 

6.  The  single  particle  albedo. 

7.  The  particle  size. 


1.1  BACKGROUND 


1.1.1  Light  Scattering  Theory 

The  scattering  of  electromagnetic  radiation  by  a  spherical  dielectric  particle  is  described  by  Mie 
theory  1  in  terms  of  the  complex  refractive  index,  m  =  m-|  -  irri2,  the  dimensionless  light-scattering  size, 

a  =  2jta /X,  and  the  angular  positionsof  the  detector,  0  and  0,  where  0  and  0  ,  are  measured  as 
indicated  in  Figure  1.  In  the  experiments  described  herein  the  detector  is  rotated  only  in  the  x-y 
plane,  so  0  =  n/2,  and  0  will  not  appear  explicitly. 


z 


Figure  1 .  The  light-scattering  coordinates. 


The  scattered  intensity,  is,  ,s  given  oy  the  following  equation 

00  +  2 

I3  =  |Z  [annn(co30)  +  bnXn(C03e)]|  ,  (1) 

n- 1 

where  the  scattering  coefficients,  an  and  bn  are  given  by 

i^(maO!/rn(cO  -  mi^n(ma)^([i) 

°n  ^(ma)^n(a)  -  mi^n(ma)tn(a)  ^ 

and  ,  t  v 

rru/-n(ma)^n(ci)  -  tfn(ma0^n(cO  ^ 

b_  ~  - .  T) 


Mie.  G  Ann.  Phys.  (Leipzig)  25.  377  (1908). 


■.  •  *  -  >  * 


The  function,  yn,  is  defined  in  terms  of  the  spherical  Bessel  function,  jn,  by  the  equation 

yn(a)  =  ctjn(a),  (' 


and  the  function,  £n,  is  given  in  terms  of  yn  and  the  Neuman  function,  Nn(a),  as  follows 

^(a)  =  Vn(a)  +  iNn(a).  (5) 

Primes  indicate  differentiation  of  the  functions  with  respect  to  their  appropriate  arguments.  The 
functions  nn(cos0)  and  Tn(cos9)  are  written  in  terms  of  Legendre  functions,  P^1),  as  follows 


n^cosG)  =  Pnn)(cos9)/sin0,  (6) 

and 

tn(C036)  =  j^[Pn0)(co3e)].  (7) 

In  general,  the  scattering  efficiency  factor,  Qsca,  depends  on  the  orientation  of  the  particle  and 
on  the  state  of  polarization  of  the  incident  light,  but  for  a  spherical  particle  Qsca  is  independent  of  both 


and  is  given  by 


_  2 


a2  n-1 


Z  (2n+ 1  )(|a n|2  +  |bn|Z>. 


The  total  extinction  efficiency  factor,  Qext,  which  consists  of  the  efficiency  of  scattering  plus 
the  absorption  efficiency,  can  be  expressed  in  terms  of  the  Mie  scattering  coefficients  as 


Qext  “  2  2  <2"+1>Re<an  +  bn>- 

n  =  1 

where  Re  refers  to  the  real  part  of  the  argument. 

Finally,  the  albedo,  A,  of  the  single  particle  is  given  by 


A  =  Qexf 


The  albedo  is  unity  for  a  particle  which  absorbs  none  of  the  electromagnetic  energy,  and  it  is  zero  for  a 
totally  absorbing  sphere. 


rr 


For  a  nonspherical  particle  the  equations  governing  the  light  scattering  are  not  tractable,  but 
many  of  the  dust  particles  encountered  in  this  research  had  reached  temperatures  high  enough  to 
melt  and,  thus,  were  spherical.  In  this  case  the  particle  size  and  its  complex  refractive  index  were 

obtained  by  comparing  theoretical  phase  functions  (ls  as  a  function  of  9)  computed  using  Equation  (1) 
with  experimental  data  for  the  scattering  intensity  as  a  function  of  angle  in  the  horizontal  plane. 


To  obtain  the  light-scattering  data  a  charged  dust  particle  was  suspended  in  a  vertically 
polarized  laser  beam  by  means  of  superposed  ac  and  dc  electrical  fields.  The  combination 
light-scattering  photometer  and  electrodynamic  balance  developed  by  Davis  and  Ray  2  was  used  for 
the  most  of  the  experiments,  but  some  experiments  were  performed  with  a  new  device  of  improved 
design.  The  former  apparatus  and  a  schematic  diagram  of  the  circuitry  are  shown  in  Figure  2. 
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Figure  2.  The  electrodynamic  balance  and  associated  circuitry. 

An  exploded  view  of  the  light-scattering  chamber  is  shown  in  Figure  3.  The  upper  and  lower 
endcap  electrodes  and  the  ring  electrode  have  inner  surfaces  of  hyperbolic  shape,  and  the  ring 
electrode  has  a  window  inserted  at  its  midplane  which  covers  the  angle  range  5’  <9<  175'.  A 
traversing  photomultiplier  tube  (PMT),  driven  by  a  stepping  motor,  was  used  to  measure  the  intensity 
of  the  scattered  light  as  it  was  rotated  around  the  light-scattering  cell.  The  output  of  the  PMT  and  the 
voltage  applied  to  the  stepping  motor  were  recorded  on  an  X-Y  plotter  and/or  stored  on  floppy  disks 
using  a  MINC-23  digital  computer  (Digital  Equipment  Corp.)  to  provide  the  record  of  intensity  versus 
angle  needed  for  comparison  with  Mie  theory. 

2  Davis,  E.  J.,  ana  Ray,  A.K.,  J.  Colloid  Surface  Sci.  75,  566  (1980). 
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Figure  3.  An  exploded  view  of  the  electrodynamic  balance. 

The  function  of  the  ac  field  generated  by  the  ring  electrode  is  to  focus  the  charged  particle  near 
the  center  of  the  balance  chamber,  and  the  dc  field  is  used  to  balance  the  gravitational  force.  In  the 
absence  of  a  dc  field  the  particle  oscillates  near  the  center  of  the  balance,  and  when  the  vertical 
component  (z-component)  of  the  dc  field  satisfies  the  following  force  balance  the  particle  remains 
stationary  at  the  null  point  of  the  chamber: 

qEdc,z  -  =  0-  (11) 

where  q  is  the  charge  on  the  particle,  E<jC  Z  is  the  vertical  component  of  the  dc  field,  m  is  the  particle 
mass,  and  g  is  the  gravitational  acceleration  constant.  The  dc  field  is  given  by 

Edc,z  =  ^o^o^zO'  (12) 


where  V0  is  the  voltage  applied  to  the  lower  endcap  electrode,  C0  is  a  geometrical  constant  of  the 
balance,  and  2z0  is  the  minimum  distance  between  the  upper  and  lower  electrodes. 

Thus,  in  the  aosence  or  other  forces  (sucn  as  aerooynamic  orag  and  photophoretic  forces)  the 
voltage,  V0,  required  to  suspend  the  particle  in  the  laser  beam  passing  through  the  center  of  the 
balance  is 


If  a  uniform  flow  of  gas  is  passed  upward  through  the  balance  chamber,  the  suspended  particle 
will  experience  a  drag  force  which  depends  on  the  gas  velocity  and  the  particle  shape.  Dahneke  3 
reviewed  the  literature  on  the  aerodynamic  drag  on  small  (low  Reynolds  numbers)  nonspherical 
bodies,  and  for  creeping  flow  shape  and  slip  correction  factors  can  be  introduced  to  give  the 
generalized  equation, 

Fd  =  K0Lcp.U00/(1+f),  (14) 

where  Fd  is  the  aerodynamic  drag  force,  K0  is  the  resistance  constant,  p.  is  the  gas  viscosity,  l_c  is  the 
characteristic  dimension  of  the  particle,  UM  is  the  gas  velocity  in  the  mainstream  of  the  flow,  and  f  is  the 
slip  correction  factor.  The  slip  correction  factor  allows  for  deviations  from  continuum  theory  which 
occur  when  the  gas  pressure  is  low.  At  normal  atmospheric  pressure  f  may  be  taken  to  be  zero,  and 
we  shall  assume  it  to  be  zero  in  this  report. 

For  spherical  particles  in  the  continuum  regime  Equation  (14)  reduces  to  Stokes’  equation, 

Fd  =  SrtapUoo,  (15) 

where  a  is  the  sphere  radius.  Equation  (15)  is  valid  for  Re  <  1,  where  the  Reynolds  number  is 
defined  by  Re  =  2aU^,'u,  and  x>  is  the  kinematic  viscosity  of  the  gas.  For  somewhat  higher  Reynolds 
numbers  Oseen’s  solution  4  can  be  used,  which  is 

Fd  =  BnaplUl  +  3Re/16  +  19Re2/2560  + . ).  (16) 

Dahneke  reported  shape  correction  factors  for  various  bodies  with  well-defined  shapes  such  as 
rods,  disks,  oblate  spheroids  and  prolate  spheroids,  but  for  aerosol  particles  of  irregular  shape 
experimental  data  are  generally  unavailable. 

To  obtain  aerodynamic  drag  data  for  the  blow-off  dust  particles  the  electrodynamic  balance  was 
modified  to  provide  a  flow  of  metered  nitrogen  upward  along  the  vertical  axis  of  the  chamber.  This 
was  accomplished  by  drilling  a  small  hole  (1 .5  mm  diameter)  in  the  bottom  electrode  to  permit  a  laminar 
jet  of  gas  to  enter  the  chamber.  The  aerodynamic  drag  force,  then,  acted  in  the  upward  direction  in 
opposition  to  gravity.  Thus,  the  dc  electrical  field  needed  to  balance  the  particle  was  altered  by  the 
flow  and  can  be  used  to  determine  the  drag  force.  When  the  jet  velocity  is  just  large  enough  to  pro- 
ouce  a  drag  torce  equal  to  the  gravitational  force,  no  dc  tieid  is  required  to  balance  the  particle,  ano  tor 

3  Dahneke,  B.  E.,  J.  Aerosol  Sci.  4,  139  (1973). 

4  Oseen,  C.,  Neure  Methcden  und  Ergebrvsse  In  der  Hydrohynamik,  Akademische  Verlag, 
Leipzig,  German.  1927 


low  Reynolds  numbers  the  critical  velocity,  Uc,  corresponding  to  this  force  balance  must  satisfy  the 
equation, 

K0LcpUc  =  mg.  (17) 


In  general,  the  force  balance  on  the  particle  is  given  by 


OQqVftJZzQ  +  Fd  -  mg  =  0.  (18) 

It  is  convenient  to  write  Equation  (18)  in  an  alternate  form  by  dividing  the  equation  by  mg  and  using 
Equation  (13)  to  eliminate  the  charge-to-mass  ratio,  q/m.  The  result  can  be  written  in  the  form 

Fd/mg  =  1  -  Vdc/V0>  (19) 

where  Vdc  and  V0  are  the  dc  voltages  required  to  suspend  the  particle  with  and  without  gas  flow, 
respectively.  Since  Vdc  is  a  function  of  the  gas  flow  rate,  the  drag  characteristics  of  a  particle  can  be 
determined  over  a  fairly  wide  range  of  flow  rates  by  simply  measuring  the  voltage  required  to  suspend 
the  particle  with  no  flow  through  the  chamber,  then  with  various  flow  rates  up  to  the  limit  of  the  stability 
of  the  particle  in  the  balance.  We  shall  report  the  drag  force  to  weight  ratio,  Fd/mg,  for  the 
experiments  performed  in  this  research. 

1 .1 .4  Calibration  of  the  Electrodvnamic  Balance 

Two  calibrations  were  involved  with  the  apparatus.  The  first  was  the  determination  of  the 
geometrical  constant,  C0,  appearing  in  the  force  balance  equations.  The  second  calibration  was  that 
of  the  jet  velocity  as  a  function  of  the  metered  gas  flow  rate. 

The  value  of  C0  was  determined  by  the  method  outlined  by  Davis  5  based  on  the  stability 
characteristics  of  the  electrodynamic  balance.  The  ac  field  which  provides  the  focusing  of  the  particle 
causes  the  particle  to  oscillate  when  the  dc  field  does  not  balance  gravity,  and  the  oscillation  is 
damped  by  the  viscous  drag  of  the  surrounding  gas.  This  is  true  provided  that  the  ac  field  strength  is 
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not  too  large  to  overcome  viscous  damping.  At  lower  frequencies  and  higher  amplitudes  of  the  ac 
voltage  the  particle  becomes  unstable  and  cannot  be  maintained  in  a  stationary  state  at  the  center  of 
the  balance.  By  measuring  the  ac  voltage  above  which  the  particle  cannot  be  stably  balanced  (the 
marginal  stability  state)  and  by  using  the  theoretical  results  of  the  stability  theory  of  the  balance,  the 
geometrical  constant  was  found  to  be  given  by  C0  =  0.471  with  a  standard  deviation  of  0.054  and  a 

mean  deviation  of  9.1%  based  on  20  experiments  performed  with  dibutyl  phthalate  droplets  of  various 
diameters. 

The  droplet  diameter  was  determined  by  light-scattering  measurements,  and  the  droplet 
mass  calculated  from  the  droplet  diameter  and  the  known  density  of  the  liquid.  The  experiments  were 
performed  by  capturing  a  droplet  in  the  balance  using  an  ac  field  lower  than  the  marginal  stability  state. 
The  dc  voltage  required  for  suspension  was  recorded,  and  light-scattering  data  were  obtained.  Then 
the  ac  voltage  was  increased  until  sudden,  large-amplitude  oscillations  were  observed  in  the 
microscope  attached  to  the  ring  electrode.  By  comparison  of  the  measured  marginal  stability  state 
with  the  theoretical  predictions  of  Davis,  the  geometrical  constant  was  determined. 

The  jet  velocity  was  calibrated  by  suspending  spherical  particles  and  droplets,  whose  sizes 
were  measured  by  light-scattering,  at  various  flow  rates.  The  critical  volumetric  flow  rate,  Qc,  required 
to  balance  gravity  and  the  no-flow  voltage,  V0,  were  measured  for  each  droplet.  The  critical  velocity, 
Uc,  corresponding  to  that  flew  rate  was  calculated  using  the  equation, 

6rtapUc(1  +  3Rec/16  +  19Rec2/2560)  -  mg  =  0,  (20) 

where  the  droplet  mass  is  obtained  from  the  droplet  radius  by  m  =  47ia2p|_/3,  where  pj_  is  the  density 
of  the  liquid.  For  Rec  <  1  this  reduces  to 

Uc  =  2pLa2g/9p.  (21) 

Equations  (20)  and  (21)  are  also  the  basis  for  determining  the  aerodynamic  radius  from  measured 
values  of  the  critical  flow  rate,  for  once  the  jet  velocity  was  calibrated  the  same  procedure  was  used  to 
determine  the  radius. 

Figure  4  shows  the  calibration  results  for  nitrogen  at  room  temperature.  At  higher  flow  rates  the 
jet  velocity  is  approximately  proportional  to  Q2,  which  is  consistent  with  laminar  jet  theory  °.  For  a  cir- 
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cular  laminar  jet,  if  the  momentum  of  the  jet  is  conserved,  we  may  write 


,oo 


J  =  2np!  0  u2rdr  =  constant, 


(22) 


where  J  is  the  flux  of  momentum,  p  is  the  fluid  density,  u  is  the  local  jet  velocity,  and  r  is  the  radial 
position.  The  velocity  distribution  for  the  laminar  jet  is  given  by 
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u  = 


_3_  _K  1 

8rf  vz  (1  +  £2/4  )2 


(23) 


where  K  is  the  kinetic  momentum  of  the  jet  (K  =  J/p),  z  is  the  distance  from  the  orifice,  and  E,  is  defined 
by 


(24) 


Figure  4.  The  calibration  of  the  jet  velocity  for  drag  measurement. 


Now  the  kinetic  momentum  can  be  approximated  by  the  momentum  of  the  gas  leaving  the  hole  in  the 


bottom  electrode.  If  the  flow  in  that  tube  is  Poiseuille  flow,  its  kinematic  momentum  is  given  by 


K  =  16Q2/37id02,  (25) 

where  Q  is  the  volumetric  flow  rate  of  gas  and  d0  is  the  diameter  of  the  hole.  Thus,  using  Equation 
(25)  in  Equation  (23)  and  evaluating  the  velocity  at  the  center  of  the  balance  (r  =  0  and  z  =  z0),  we 
obtain 


u0  =  2Q2/7t2d02\)Z0. 


(26) 


Equation  (26)  is  also  plotted  on  Figure  4.  The  experimental  data  for  high  flow  rates  fall  below 
but  are  parallel  the  theoretical  velocities  predicted  using  Equation  (26).  At  low  flow  rates  the  velocities 
diverge  significantly  from  the  results  predicted  for  a  free  laminar  jet,  for  the  jet  velocities  are  greater 
than  predicted  values.  This  is  a  result  of  jet  expansion  and  entrainment  of  surrounding  gas  by  the  jet, 
as  indicated  in  Figure  5a  and  5b.  Figure  5a  shows  the  jet  velocity  profiles  for  two  vertical  planes 
predicted  using  Equation  (23)  for  the  free  jet  for  a  nitrogen  flow  rate  of  10  ml/min,  and  Figure  5b 
shows  profiles  for  a  flow  rate  of  30  ml/min.  At  the  lower  flow  rate  the  jet  expands  to  interact  with  the 
walls  of  the  container  before  it  reaches  the  midplane,  so  the  je ‘  cannot  be  considered  to  be  a  free  jet. 
At  the  higher  flow  rate  there  is  less  interaction  with  the  bounding  surfaces,  and  the  jet  approaches 
free-jet  behavior. 
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Figure  5  a.b.  Calculated  jet  velocity  profiles  for  a  nitrogen  flow  rate  of  10  ml/min  and  30  ml/mi, 
respectively 


SECTION  2 

EXPERIMENTAL  PROCEDURES  AND  DATA  ANALYSIS 


2.1  EXPERIMENTAL  TECHNIQUES 

2.1.1  Particle  Injection 

Liquid  droplets  and  the  blow-off  dust  particles  were  charged  and  injected  into  the  balance  in 
the  same  manner.  A  small  liquid  drop  was  partly  formed  at  the  tip  of  a  stainless  steel  capillary,  and  a  dc 
voltage  of  2  to  5  kV  was  suddenly  applied  to  the  capillary  from  a  high  voltage  supply.  The  high  charge 
density  on  the  liquid  surface  reduced  the  effective  surface  tension  to  zero,  causing  the  droplet  to  fall 
into  the  balance.  With  the  ac  voltage  preset  to  approximately  1000  V  rms,  the  droplet  was  trapped. 
The  dc  voltage  was  then  adjusted  to  stop  the  droplet  oscillation,  as  observed  through  a  70X 
microscope  attached  to  the  balance.  The  droplet  was  suspended  at  the  null  point  through  which 
passed  a  vertically  polarized  laser  beam.  For  most  of  the  experiments  a  tunable  dye  laser  (Coherent 
Inc.  Model  CR-599),  pumped  by  an  argon-ion  laser  (Coherent  Inc.  Model  Innova  90-5),  was  the  light 
source,  but  in  November  1984  the  pump  laser  failed,  and  the  laser  system  was  replaced  by  an 
available  Spectra  Physics  laser  operating  at  a  fixed  wavelength  of  632.8  nm.  Rhodamine  6G  was  the 
dye  used  in  the  dye  laser,  and  this  provided  wavelengths  in  the  range  575-650  nm.  Most  of  the 
light-scattering  data  reported  herein  were  obtained  at  605  nm,  and  data  in  the  report  on  Phase  II  of  the 
program  were  for  wavelengths  of  457.6.  488.0  and  514  5  nm.  Due  to  failure  of  the  pump  laser  it  was 
not  possible  to  make  measurements  in  the  near-ultra  violet  region  of  the  spectrum. 

When  blow-off  particles  were  studied,  the  particles  were  washed  from  the  filters  or  wiping 
material  supplied  by  SAI,  and  a  droplet  consisting  of  a  dilute  suspension  of  the  particles  was  injected 
into  the  balance  The  water  evaporated  within  a  few  seconds,  leaving  the  charged  particle  in  the  laser 
beam. 

2.1.2  Light-Scattering  Data  Acquisition 

Prior  to  capturing  a  particle  the  light-scattering  data  acquisition  system  was  activated  to  record 
background  intensities.  Background  signals  were  also  recorded  at  the  end  of  a  run.  These  data  were 
used  to  correct  the  ugnt-scattering  data  obtained  with  a  suspended  particle. 

After  a  particle  had  been  suspended  in  a  dry  gas  stream  for  about  two  hours  to  remove  all  water 
from  its  surface  the  rotating  photomultiplier  tube  was  activated  and  the  scattering  intensities  versus 


angle  were  recorded  on  an  X-Y  plotter.  Based  on  the  plotted  data,  adjustments  to  the  detector 
voltage  or  other  desired  corrections  were  made,  and  then  the  data  were  digitized  and  stored  on  floppy 
disks  of  a  digital  computer.  The  stored  data  were  later  replotted  using  a  Hewlett-Packard 
printer/plotter  (Model  7245B)  connected  to  the  computer. 

Usually  a  forward  scan  (increasing  ,)  and  a  backward  scan  (decreasing ,)  with  the  PMT  were 
made  and  compared.  If  a  particle  was  nonspherical,  it  was  usually  observed  to  "twinkle"  in  the 
microscope,  and  the  forward  and  backward  scans  would  not  match  because  of  particle  rotation. 

A  representative  experimental  phase  function  for  a  blow-off  dust  particle  is  shown  in  Figure  6. 
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Figure  6.  A  typical  phase  function  obtained  for  a  blow-off  dust  particle. 


After  the  light-scattering  data  were  ol  ned,  the  gas  (low  rate  was  varied  and  the  dc  voltage 
was  recorded  as  a  function  of  flow  rate.  In  many  runs  the  tlow  rate  was  increased  until  the  particle  was 
blown  from  the  balance.  At  higher  flow  rates  a  particle  would  begin  to  move  off  the  centerline  of  the 
balance,  and  at  yet  higher  flow  rates  it  would  be  lost 


Typical  data  for  dust  particles  are  presented  in  Figure  7  as  a  plot  of  the  dc  voltage  required  for 
stable  suspension  versus  the  nitrogen  flow  rate.  In  this  case  the  particles  were  positively  charged, 
and  the  voltage  decreased  nonlinearly  from  some  positive  value,  crossed  the  axis  at  the  critical  flow 
rate,  Qc,  and  continued  to  decrease  as  the  flow  rate  was  increased. 
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Figure  7.  The  dc  suspension  voltage  as  a  function  of  nitrogen  flow  rate  for  representative  blow-off 
dust  panicles. 

The  data  from  Figure  7  are  replotted  in  terms  of  'he  drag  force-to-weight  ratio.  Fw/mg,  versus 
nitrogen  velocity  in  Figure  8.  For  each  particle  the  drag  force  is  seen  to  be  a  linear  function  of  velocity, 
which  indicates  that  the  flow  was  in  the  creeping  flow  regime.  Note  that  for  Sample  136H  the  ratio 
F^/mg  exceeded  1 1 0,  and  yet  the  drag  force  varied  linearly  with  the  jet  velocity. 
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Figure  8.  The  drag  force-to-weight  ratio  as  a  function  of  nitrogen  velocity  for  the  data  of  Figure  7. 


An  atypical  set  of  data  for  three  particles  is  shown  in  Figures  9  and  10.  These  particles  were  of 
irregular  shape  and  show  nonlinear  drag  force  behavior.  It  appears  that  particle  60M  underwent 
realignment  as  the  gas  velocity  increased,  for  it  is  clear  from  Figure  10  that  the  slope  of  the  curve  of 
Fjj/mg  versus  U  changed  significantly  as  the  velocity  increased.  Sample  60P  and,  to  a  lesser  extent, 
Sample  69K  also  showed  systematic  deviations  from  linearity.  Particle  69K  had  an  estimated 
aerodynamic  radius  of  only  0.582  pm,  which  was  one  of  the  smaller  particles  examined,  and  the  drag 
force-weight  ratio  reached  900  before  the  particle  was  lost. 
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Figure  1 0.  Nonlinear  drag  characteristics  of  blow-off  dust  particles 


2.2  INTERPRETATION  OF  THE  LIGHT-SCATTERING  DATA 

The  analysis  of  the  light-scattering  phase  functions  follows  that  of  Davis  and  Ravindran  7  who 
applied  the  Marquardt-Levenberg  algorithm  ®  to  determine  the  complex  refractive  index  and  optical 
size  from  data  by  comparison  with  Mie  theory. 

Forty  equally-spaced  data  points  in  the  angle  range  30*  to  130*  were  selected,  and  the 
intensities,  lt,  predicted  from  Mie  theory  for  these  angles  were  computed  using  trial  values  of  the 

unknown  parameters  (m  =  m-j-im2,  and  a).  The  voltages  from  the  photodetector  circuit  were 
converted  to  normalized  intensities  using  the  relation, 

le(9i)  =  oIV(0j)  -  Vb(eij)],  (27) 


where  cr  is  a  scale  factor,  V(0j)  is  the  signal  at  angle  9  j,  V£,(0j)  is  the  background  signal  for  that  angle, 
and  le(0j)  is  the  experimental  intensity  at  0j. 

The  scale  factor  was  considered  to  be  one  of  the  parameters  in  the  optimization  procedure. 
The  procedure  was  to  minimize  the  function, 


F(a.m,o)-SjWi|Rj|2 


(28) 


where  Rj  is  the  residual  defined  by 


Ri  =  fe<0i>  -  ft(0i>- 


(29) 


and  Wj  is  a  weighting  function,  which  was  taken  to  be  Wj  =  sin20j. 

To  minimize  the  objective  function,  F(a,m,o),  the  Marquardt-Levenberg  routine  performs  an 
optimum  interpolation  between  the  Taylor  series  and  gradient  (steepest  descent)  methods.  The 
computer  code  was  run  with  various  trial  values  in  an  attempt  to  find  the  global  minimum  of  the 
function.  To  speed  up  the  calculations  the  trial  value  of  the  optical  size,  a,  was  estimated  from  the 
measured  aerodynamic  size,  which  was  obtained  from  the  drag  force  experiments.  A  good  trial  value 
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of  owas  obtained  by  rescaling  the  experimental  maximum  peak  and  minimum  trough  values  to  the 
maximum  and  minimum  values  for  the  angle  range  computed  from  Mie  theory.  Usually,  m2,  the 
imaginary  component  of  the  refractive  index,  was  set  equal  to  zero  to  start  the  computations. 

Once  the  optimum  values  of  a,  m  and  a  were  computed  the  phase  function  was  calculated  from' 
Mie  theory  and  compared  with  the  experimental  data  to  examine  the  quality  of  the  "fit".  We  shall 
present  the  results  below. 

An  alternate,  and  simpler  procedure,  for  determining  a  trial  value  of  the  light-scattering  size  was 
developed,  which  yielded  accurate  sizes  for  larger  particles  (a  >  2  pm)  and  a  reasonably  good  estimate 
of  size  for  smaller  particles.  This  technique  was  based  on  Mie  theory  computations  for  spheres  of 
various  sizes  and  refractive  indices,  which  indicate  that  for  large  a  the  phase  functions  are  relatively 
insensitive  to  refractive  index  but  show  a  monotonically  increasing  number  of  peaks  for  a  given  angle 
range,  say  from  40*  to  90*.  Figure  1 1  presents  the  results  of  phase  function  computations  for  several 
refractive  indices  in  the  range  of  interest  plotted  as  the  number  of  peaks  per  50*  versus  the  optical 
size,  a.  Thus,  simply  by  counting  the  number  of  peaks  over  a  50*  range  it  was  possible  to  determine 
a  within  about  3%  for  a  >  20. 
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Figure  1 1 .  The  number  of  phase  function  peaks  per  50*  as  a  function  of  a. 


SECTION  3 
RESULTS 


3.1  SAMPLE  IDENTIFICATION 

The  identification  of  the  blow-off  dust  samples  follows  the  notation  used  by  SAI  for  data 
collection  in  the  solar  furnace  experiments.  Table  1  lists  the  run  numbers,  soil  type  and  a  description 
which  indicates  where  the  sample  was  taken  from,  for  example,  from  an  aerosol  filter  or  from  wall 
wipings.  We  shall  identify  each  of  our  experiments  by  using  the  SAI  run  number  together  with  an 
alphabetic  character  starting  with  A.  For  example,  Sample  136C  represents  the  third  dust  particle 
examined  from  SAI's  run  136. 

TABLE  1 .  Identification  of  samples. 


SAI  Run  Number 

Soil  Type  Code 

Description 

57 

Silt  P-14 

Wall  Wipings 

60 

Clay  P-13 

Wall  Wipings 

69 

RB5I-2 

Wall  Collection 

104 

RC5Y-2 

Wall  Wipings 

107 

RC4U-2 

Wall  Wipings 

136 

Soil  RBU2#15 

Wall  Wipings 

138 

Black  Tile 

Wall  Wipings 

1180 

Wall-5  F 

Wall  Wbings 

3.2  PARTICLE  CHARACTERIZATION 

Each  macroscopic  sample  of  blow-off  dust  was  examined  in  the  scanning  electron  microscope 
to  determine  qualitatively  the  size  distribution  and  shapes  of  particles.  The  microscope  was  equipped 
with  energy  dispersion  x-ray  spectrometry  (EDAX)  capabilities  for  elemental  analysis.  Figures  12a,b 
through  19a,b  show  repesentative  results  from  the  scanning  electron  microscope  studies.  We  note 
that  in  most  of  the  samples,  with  the  possible  exception  of  Runs  136  and  138,  spherical  particles  are 
observed,  indicating  that  the  temperature  in  the  solar  furnace  was  high  enough  to  melt  many  of  these 
mixed  silicates. 

The  EDAX  analyses  indicate  that  the  samples  were  mostly  mixed  silicates  of  aluminum,  calcium, 
magnesium  and  potassium.  Iron  was  found  in  significant  amounts  in  samples  from  Runs  60  and  136. 
Table  2  lists  the  results  of  EDAX  analysis. 
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TABLE  2.  Results  of  EDAX  analysis. 


Run  Number 


r 

r 

r 

i 


i 


57 

60 

69 

104 

107 

136 

138 

1180 


Elements  Present  in  Sample 

Al,  Ca,  K,  Na,  Si,  Fe  (trace) 

Al,  Ca,  Fe,  K.  Si,  Ag  (trace) 

Al,  Ca,  K,  Na,  Si,  traces  of  Ag,  Fe,  Cu 
Al,  Si,  Cu  (trace) 

Al,  Ca.  Si,  traces  of  Cu,  Fe 
Al,  Ca,  Fe,  K,  Cu  (trace) 

Ag,  Al,  Ca,  Si,  Fe  (trace) 

Al,  Ca,  Si,  Cu  (trace) 


It  is  not  possible  to  determine  the  types  of  minerals  and  their  chemical  formulas  in  the  samples 
from  the  EDAX  analyses  because  the  EDAX  results  are  qualitative,  but  it  is  reasonable  to  assume  that 
minerals  such  as  Andesine  [(Ca0,Na20)Al2C>3'4Si023,  Anorthite  (Ca0Al203-2Si02),  Anorihoclase 
[(Na,K)20-Al203'6Si02)].  Laboradorite  (NaAISisOg  +  CaA^S^C^),  Orthoclase 
(K20-Al203'6SiC>2),  Kaolinite  (A^C^^SiC^F^O),  Gehlenite  (Ca2Al2SiC>7)  or  other  mixed  silicates 
were  present.  The  measured  refractive  indices  are  consistent  with  the  bulk  refractive  indices  of  such 
minerals,  as  we  shall  show. 


3.3  AERODYNAMIC  DRAG  CHARACTERISTICS 


The  ratio  of  the  aerodynamic  drag  force,  Fd,  to  the  weight,  mg,  of  a  particle  was  obtained  from 
the  dc  voltage  data  using  Equation  (19).  The  aerodynamic  radius  was  calculated  from  the  critical 
velocity,  Uc  required  for  the  balance  of  the  drag  force  and  the  gravitational  force  using  Equation  (20) 

or  (21).  Figures  20-27  show  the  ratio  Fd/mg  versus  particle  Reynolds  number  for  numerous  dust 
particles  of  various  sizes  and  shapes  from  various  soil  samples.  The  particle  Reynolds  numbers  were 
calculated  using  the  measured  aerodynamic  radii  Because  of  the  small  sizes  of  the  particles  (a  <  20 
pm)  the  Reynolds  numbers  did  not  exceed  0.25,  so  all  of  the  data  were  in  the  creeping  flow  regime. 
Figures  20-27  show  that  the  drag  forces  are  linear  with  respect  to  the  Reynolds  numbers.  Additional 
results  are  provided  in  the  Appendix. 

To  examine  the  effects  of  particle  shape  and  roughness  on  the  aerodynamic  drag  it  is 
convenient  to  report  the  drag  force  in  terms  of  the  particle  mobility,  B,  defined  by 


B  =  lUFd. 


(30) 


For  a  spherical  particle  in  the  Stokes  regime  of  creeping  flow  the  mobility  is  given  by 

B  =  1/6?rap,  (31) 

but  for  nonspherical  particles  B  depends  on  the  shape  and  roughness. 

The  drag  force  can  be  calculated  from  the  ratio  F^/mg  as  a  function  of  gas  velocity  provided  that 
we  can  determine  the  particle  mass.  The  slope  of  a  plot  of  F^/mg  versus  U  is  given  by  l/Bmg.  Thus, 
the  reciprocal  of  the  mobility  is  obtained  by  multiplying  the  measured  slope  by  the  particle  weight,  mg. 
For  Figures  20-27,  which  have  as  the  abscissae  the  Reynolds  number  rather  than  the  velocity  U,  the 
slopes  of  the  lines  are  inversely  proportional  to  the  particle  mobility. 

To  calculate  the  mobility  from  the  data  on  Fd/mg  versus  U  the  particle  mass  was  estimated  by 

assuming  the  particle  to  have  the  volume  of  a  sphere  having  the  same  aerodynamic  radius  and  a 
density  equal  to  2.855  g/cm3,  the  average  density  of  silicates  with  the  elemental  composition  and 
refractive  indices  of  the  particles  examined  in  this  study.  Table  3  summarizes  the  mobilities  calculated 
in  this  way  compared  with  the  Stokesian  mobilities  obtained  using  the  aerodynamic  radii  substituted  in 
Equation  (31).  For  spherical  particles  the  two  mobilities  should  be  identical  if  the  assumed  density  is 
correct. 
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For  most  of  the  data  listed  in  Table  3  the  mobilities  based  on  Stokesian  drag  do  not  differ 
greatly  from  the  mobilities  computed  from  the  drag  force  data,  but  there  are  a  few  exceptions. 
Particles  57A,  57F,  60J,  1180B  and  1180F  had  mobilities  at  least  30%  lower  than  those  calculated 
from  Stokes'  law,  and  several  other  particles  had  mobilities  at  least  20%  lower.  Some  of  these 
discrepancies  are  due  to  uncertainties  in  the  particle  size  and  density. 

Many  of  the  particles  listed  as  irregularly  shaped  in  Table  3  had  mobilities  close  to  those 
predicted  using  Stokes'  law,  for  example,  irregular  particles  60G,  60K,  60T,  60U,  69C,  69E,  69G,  691, 
1048,  107E,  136B,  136C,  136F,  138E,  1180C,  1180D  and  1180E.  From  this  result  it  is  concluded 
that  most  of  the  particles  were  sufficiently  similar  to  spheres  to  have  the  drag  characteristics  of 
spheres. 
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Figure  23.  The  ratio  of  the  arag  torce  to  particle  weignt  as  a  function  of 
the  particle  Reynolds  number  for  particles  from  Sample  104 
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Figure  26.  The  ratio  of  the  drag  force  to  particle  weight  as  a  function  of 
the  particle  Reynolds  number  for  particles  from  Sample  138. 
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Figure  27.  The  ratio  of  the  drag  force  to  particle  weight  as  a  function  of 

the  particle  Reynolds  number  for  particles  from  Sample  1180. 


TABLE  3.  The  radii  and  mobilities  of  dust  particle. 


Sample 

Aerodynamic 

Light  Scattering 

Stokes  Mobility, 

Mobility, 

Shape 

Radius, 

Radius, 

Sq.  (31) 

Eq. (30) 

um 

pm 

aftxlQ'6 

s/QX  IQ'6 

57A 

3.78 

3.32 

0.915 

0.634 

Spherical 

57B 

4.46 

- 

0.303 

0.274 

Irregular 

57C 

3.31 

3.52 

0.920 

0.904 

Spherical(?) 

57D 

4.30 

3.91 

0.707 

0.715 

Spherical 

57E 

3.39 

3.56 

0.897 

0.873 

Spherical(?) 

57F 

2.67 

2.44 

1.14 

0.689 

Spherical 

57G 

10.04 

9.87 

0.303 

0.274 

Spherical 

60A 

- 

4.24 

0.718 

- 

Spherical 

60B 

1.56 

1.79 

1.95 

1.960 

Spherical 

60C 

3.51 

3.76 

0.867 

0.610 

Irregular 

60D 

1.67 

2.21 

1.82 

1.83 

Spherical 

60E 

3.33 

2.60 

0.914 

0.891 

Spherical 

60F 

3.18 

- 

0.957 

0.886 

Irregular 

60G 

2.72 

- 

1.116 

1.12 

Irregular 

60H 

4.18 

- 

0.726 

0.690 

Irregular 

601 

4.30 

3.07 

0.707 

0.786 

Spherical(?) 

60J 

4.10 

2.59 

0.7s  . 

0.473 

Irregular 

60K 

1 .66 

- 

1.83 

1.71 

Irregular 

60L 

1.12 

2.17 

2.72 

2.54 

Spherical 

60M 

1.79 

- 

1.70 

- 

Irregular 

60N 

1.99 

- 

1.53 

- 

Irregular 

60P 

2.42 

- 

1.26 

- 

Irregular 

60Q 

12.6 

- 

2.43 

- 

Fiber(?) 

60R 

0.394 

- 

7.72 

8.50 

Irregular 

60S 

1.13 

- 

2.70 

- 

Irregular 

60T 

4.67 

- 

0.651 

0.650 

Irregular 

60U 

0.792 

- 

3.84 

3.61 

Irregular 

69A 

5.42 

- 

0.560 

0.668 

Irregular 

69B 

183 

- 

0  166 

0.126 

Aggregate^: 

69C 

3.17 

5.27 

0  577 

0  533 

Spherical(?) 

e!g 

Aerodynamic 

TABLE  3.  T 

Light  Scatte 

Radius, 

Radius, 

_UHL 

pm 

D 

4.35 

- 

E 

1.49 

1.48 

iF 

6.03 

6.24 

iG 

4.60 

- 

IH 

2.51 

2.60 

ll 

3.12 

- 

)J 

3.50 

3.42 

)K 

0.582 

- 

\A 

1.43 

1.55 

IB 

3.85 

- 

■C 

3.09 

2.60 

.D 

3.96 

3.76 

IE 

3.85 

3.62 

7A 

3.61 

1.33 

7B 

6.86 

1.40 

’C 

5.80 

6.00 

7D 

5.23 

1.36 

7E 

5.23 

- 

5A 

2.74 

2.32 

5B 

3.52 

- 

3C 

3.15 

- 

3D 

3.85 

2.84 

6E 

3.31 

3.47 

6F 

4.96 

- 

6G 

4.61 

6.24 

3H 

1.73 

- 

8A 

6.10 

5.78 

8B 

3.26 

3.90 

8C 

1.47 

1.32 

i  and  mobilities  of  dust  particle  (Continued) . 


Stokes  Mobility, 
Eq.pl) 


Mobility, 
Eq.  (30) 


Shape 


S/gxJir6 

0.699 

S/QX  IQ'6 

0.783 

Rodlike(?) 

2.05 

2.05 

Spherical 

0.504 

0.499 

Spherical 

0.660 

0.717 

Irregular 

1.21 

1.23 

Spherical 

0.974 

1.09 

Irregular 

0.868 

0.867 

Spherical 

0.522 

- 

Irregular 

2.13 

2.19 

Spherical 

0.789 

0.810 

Irregular 

0.984 

0.835 

SphericalP) 

0.763 

0.757 

Spnerical(?) 

0.789 

0.786 

Spherical 

0.842 

0.608 

Spherical 

0.443 

0.449 

Spherical 

0.507 

0.520 

Spherical 

0.532 

0.417 

Spnericai 

0.582 

0.551 

Irregular 

1.11 

1.13 

Spherical 

0.864 

0.778 

Irregular 

0.936 

0.852 

Irregular 

0.789 

0.847 

Spherical 

0.920 

0.781 

Spherical 

0.612 

0.612 

Irregular 

0.660 

0.608 

Spherical 

1.76 

1.35 

Irregular 

0.495 

0.584 

Spherical 

0.932 

0.775 

Spherical 

2.07 

2.08 

Spherical 

TABLE  3.  The  radii  and  mobilities  of  dust  particle  (Concluded). 


Sample 

Aerodynamic 

Light  Scattering 

Radius, 

Radius, 

urn 

urn 

138D 

1.36 

2.52 

138E 

0.428 

- 

138F 

4.08 

4.23 

138G 

6.63 

- 

138H 

14.0 

- 

1381 

1.88 

1.86 

1180A 

4.18 

3.91 

1180B 

3.02 

2.52 

1180C 

3.59 

- 

1180D 

3.43 

- 

1180E 

3.31 

- 

1180F 

3.04 

- 

Stokes  Mobility, 

Mobility, 

Shape 

Ed-  (31) 

Eq.  (30) 

s/a  x  1 0‘6 

s/a  x  1 0'6 

2.24 

2.95 

Spherical(?) 

7.10 

7.05 

Irregular 

0.745 

0.744 

Spherical 

0.458 

0.333 

Irregular 

0.217 

0.288 

Rod(?) 

1.62 

1.62 

Spherical 

0.727 

0.713 

Spherical 

1.01 

0.654 

Spherical 

0.844 

0.885 

Irregular 

0.886 

0.885 

Irregular 

0.919 

0.918 

Irregular 

0.999 

0.410 

Irregular 

RELATIVE  INTENSITY 


3.4  LIGHT  SCATTERING  RESULTS 


The  main  light  scattering  results  were  reported  under  Phase  II  of  the  research  program,  but 
additional  data  were  obtained  at  a  wavelength  of  605  nm  in  Phase  III.  The  phase  function  data  for  all 
particles  studied  were  analyzed  using  the  optimization  algorithm  discussed  in  Section  2.2,  and  typical 
results  are  plotted  in  Figures  28-34.  The  experimental  intensities  are  plotted  for  a  discrete  set  of  data 
points,  and  the  "best  fit"  results  from  Mie  theory  are  plotted  as  continuous  curves  on  the  figures. 
Additional  comparisons  are  made  in  Appendix  B. 

The  agreement  between  the  experimental  data  and  Mie  theory  is  not  as  good  as  is  usually 
obtained  for  liquid  droplets,  which  are  perfect  spheres,  for  the  dust  particles  are  not  likely  to  be 
perfectly  spherical  even  if  melted.  Furthermore,  the  particles  might  be  somewhat  inhomogeneous 
due  to  partial  crystallization,  and  that  would  affect  the  light  scattering.  In  general,  the  positions  of 
peaks  and  troughs  match  the  computed  phase  functions,  but  the  amplitudes  do  not  match  particularly 
well.  Table  4  summarizes  the  findings  for  those  particles  for  which  the  optimization  program 
converged. 


Figure  28.  A  comparison  between  light-scattering  data  for  Sample  57F 
and  Mie  theory  for  a  =  25.3  and  m  =  1 .4242  -  id. 0057. 
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RELATIVE  INTENSITY 


angle  in  degrees 


ANCLE  IN  OEGREES 


Table  4.  Light  scattering  results  for  dust  particles  at  a  wavelength  of 
605  nm  obtained  by  the  optimization  procedure. 


a  Number 

& 

Radius,  urn 

Refractive  Index 

7A 

34.5 

3.32 

1.4734 -i0.0081 

7D 

40.6 

3.91 

1.4574  -i0.0043 

7F 

25.3 

2.44 

1.4242  -  i0.0057 

>0B 

18.6 

1.79 

1.4157 -i0.0079 

>0D 

22.9 

2.21 

1.4475 -i0.0009 

JOE 

27.0 

2.60 

1.6469  -  i0.0038 

501 

31.8 

3.07 

1.5238  -  iO .001 8 

50J 

26.9 

2.59 

1.5725  -  i0.0002 

59Ea 

14.7 

1.48 

1.4832 -i0.0060 

59Ha 

25.8 

2.60 

1.4238  -  i0.010 

04A 

16.1 

1.55 

1.3908  -i0.0091 

I04E 

37.6 

3.62 

1.6175  -  i0-0015 

107A 

13.8 

1.33 

1.6044  -i0.0093 

107B 

14.5 

1.40 

1.6049  -  i0  0042 

107D 

14.1 

1.36 

1.5457-10.0042 

136A 

24.1 

2.32 

1.5154 -i0.0007 

136D 

29.5 

2.84 

1.5047  -  S0.0107 

136E 

36.1 

3.47 

1.5447 -i0.0018 

138B 

40.5 

3.90 

1.5591  -  i0.0022 

138C 

13.7 

1.32 

1.3948 -iO.0 169 

1180A 

40.6 

3.91 

1.5177 -i0.0045 

1180B 

26.2 

2.52 

1.4689  -i0.0030 

Measured  at  a  wavelength  of  632.8  nm. 


3.5  DISCUSSION  OF  RESULTS 

The  light-scattering  data  indicate  that  the  real  components  of  the  refractive  index,  m-j,  for  a 
wavelength  of  605  nm  are  in  the  range  1 .3908-1 .6469,  which  is  consistent  with  the  bulk  values  of  the 
refractive  index  of  various  mixed  silicates.  For  example,  for  the  sodium  line  (X  =  589.3  nm)  Andesine 
has  a  refractive  index  of  1.544-1.563,  Anorthite  1.577-1.590,  Anorlhoclase  1.523-1.529,  Orthoclase 
1 .518-1 .539,  Kaolinite  1 .522-1 .570,  and  Gehlenite  1 .658-1 .669.  Data  are  not  available  with  which  to 
compare  the  imaginary  component  of  the  refractive  index,  but  Mie  theory  computations  indicate  that 
the  predicted  phase  function  is  not  highly  sensitive  to  the  imaginary  component  for  the  low  range  of 
values  estimated  here.  The  main  effect  of  the  imaginary  component,  m2,  is  to  alter  the  scale  factor  or 
amplitude  of  the  phase  function  and  not  the  angular  positions  of  peaks  and  troughs.  Thus,  the 
uncertainties  in  the  reported  values  of  m2  are  great,  and  the  calculated  albedos  are  not  accurate.  It 
should  be  concluded  that  the  albedos  are  of  order  one,  and  there  is  essentially  no  significant  energy 
absorption  by  the  dust  particles  in  the  range  of  visible  light. 

The  scanning  electron  micrographs  show  that  for  Samples  57,  60,  69,  104,  107  and  1180 
many  spherical  particles  were  generated.  Even  those  particles  which  had  not  melted  to  become 
spherical  were  more  sphere-like  than  rod-like  or  disk-like,  so  their  aerodynamic  drag  characteristics 
might  be  expected  to  approximate  those  of  spheres.  The  experimental  data  ana  calculated  mobilities 
indicate,  indeed,  that  the  particles  behaved  like  spheres.  Samples  136  and  138  (see  Figures  17a  and 
18a)  showed  a  wide  variety  of  irregularly  shaped  particles,  but  the  measured  mobilities  do  not  differ 
greatly  from  the  Stokes  mobility  for  spheres. 

For  the  range  of  gas  velocities  which  could  be  used  with  the  electrodynamic  balance,  the 
Reynolds  numbers  were  generally  in  the  creeping  flow  regime,  and  the  measured  ratios  F^/mg  were 
found  to  be  linear  functions  of  Reynolds  number,  which  is  characteristic  of  creeping  flow.  Very  few 
particles  were  observed  to  have  nonlinear  drag  characteristics  (see  Figure  10),  and  those  which  were 
not  typical  were  probably  extraneous  particles  such  as  small  fibers  from  the  filters  or  wall  wiping 
material.  It  is  to  be  expected  that  at  higher  Reynolds  numbers  all  of  the  particles  would  have  a 
nonlinear  Reynolds  number  dependence. 

For  purposes  of  estimating  the  properties  of  blow-off  dust  particles  it  is  recommended  that  the 
bulk  refractive  index  be  used  for  the  types  of  minerals  expected,  that  the  single  particle  albedos  be 
considered  to  be  unity,  and  that  the  aerodynamic  drag  characteristics  be  taken  as  that  of  a  sphere  of 
the  same  volume  as  the  dust  particle  in  question. 
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Figure  40.  Drag  force  data  for  Sample  136. 
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Figure  41 .  Drag  force  data  for  Sample  138. 
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